Several Small Josephson Junctions in a Resonant Cavity: 
Deviation from the Dicke Model 
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We have studied quantum-mechanically a system of several small identical Josephson junctions 
in a lossless single-mode cavity for different initial states, under conditions such that the system 
is at resonance. This system is analogous to a collection of identical atoms in a cavity, which 
is described under appropriate conditions by the Dicke model. We find that our system can be 
well approximated by a reduced Hamiltonian consisting of two levels per junction. The reduced 
Hamiltonian is similar to the Dicke Hamiltonian, but contains an additional term resembling a 
dipole-dipole interaction between the junctions. This extra term arises when states outside the 
degenerate group are included via degenerate second-order (Lowdin) perturbation theory. As in 
the Dicke model, we find that, when N junctions are present in the cavity, the oscillation frequency 
due to the junction-cavity interaction is enhanced by y/N. The corresponding decrease in the Rabi 
oscillation period may cause it to be smaller than the decoherence time due to dissipation, making 
these oscillations observable. Finally, we find that the frequency enhancement survives even if the 
junctions differ slightly from one another, as expected in a realistic system. 

PACS numbers: 03.65.Ud, 03.67.Lx, 42.50.Fx, 74.50.+r 



I. INTRODUCTION 

In a previous paper,0 we numerically studied a small 
Josephson junction capacitively coupled to the single 
mode of a resonant cavity. We found that the quantum 
states of the combined system are strongly entangled only 
at certain special values of the gate voltage across the 
junction, corresponding to a one- or two-photon absorp- 
tion or emission process by the junction Cooper pairs. 
We also showed that at a gate voltage corresponding to 
the one-photon process, our results, obtained from exact 
diagonalization, agree well with those abtained from a 
reduced Jaynes-Cummings HamiltonianEI which includes 
only two junction states and two photon states. 

In the present paper, we extend this work to several 
similar junctions placed in a resonant cavity within a 
radiation wavelength. Each junction is assumed to be 
voltage biased by an independent voltage source, but all 
bias voltages are assumed to be equal. As in the single- 
junction case, we solve for the eigenstates of this coupled 
system by direct diagonalization, and compare the results 
to those obtained from a reduced Hamiltonian which in- 
cludes only two energy levels per junction. We also com- 
pare thfi interaction to that in a system first discussed by 
Dicke,cl consisting of an assembly of identical two-level 
atoms in a single-mode cavity. 

Before proceeding further, we note some distinctions 
between the system considered here and anotheij-|System 
which has recently been studied experimentallyfl These 
experiments involve a two dimensional array of junctions, 
such that each junction is coupled to a single-mode reso- 
nant cavity. For suitable junction and cavity parameters, 
the cavity interaction causes the array to phase-lock and 
to radiate coherentlK.|iaJ|0 the cavity. The models used to 
describe such arraysoEnZllj typically involve a semiclassical 
approximation, appropriate for large numbers of photons 



in the cavity. By contrast, our system involves only a few 
photons, and needs to be studied quantum mechanically. 
Another distinction is that, in Ref. |], the junctions are 
biased onto either the resistive or the superconducting 
branch of the I-V characteristic, whereas in the present 
model, only the latter branch is considered. 

The present work may also be relevant to quantum 
computation^ Josephson devices have been proposed as 
possible quantum bits (qubits) in computing systems, in 
which quantum logic operations could be perfarmod .by. 
controlling gate voltages or magnetic fields. B^Mii3llj'L3 
Thus far, two alternative realizations of qubits have been 
proposed for Josephson junctions: a charge qubit and a 
flux qubit. In the first case, the charging energy of the 
junction is much larger than the Josephson energy, so 
that the charge on one of the superconducting islands is 
nearly a good quantum number. Qubits of this type-have 
been studied experimentally by Nakamura et al. .E2MIII In. 
the second case, studied, experimentally by Mooij et al. E£l 
and Feigel'man et al. ,0-3 the Josephson energy predomi- 
nates, and the phase difference across the junction is well 
defined. In the present paper, the system is constructed 
from junctions in which the charging and Josephson en- 
ergies are comparable in magnitude; in other respects, it 
somewhat resembles that studied in Ref. 

The findings of this paper.-^nay also be relevant for 
studying "quantum leakage"l!3 in Josephson-junction- 
based multi-qubits. Quantum leakage occurs when the 
computational space of the physical system is a subspace 
of a larger Hilbert space. Such leaking is an intrinsic 
source of decoherence, even in the absence of dissipation, 
and imposes a time limit beyond which the system can no 
longer perform quantum computations. Thus, comparing 
the results of numerically diagonalizing the full Hamilto- 
nian and the reduced Hamiltonian provides a means of 
studying this decoherence. Such a comparison was done 
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in Ref. [15|, and is also done here for a different Hamilto- 
nian and a different parameter region. 

Finally, because of the analogy between our system and 
a group of two-level atoms in a single-mode cavity, the 
present work may be relevant to the field of cavity quan- 
tum electrodynamics.Ej Specifically, using the system 
considered here, it may be possible to "engineer" a so- 
called Einstein-Podolsky-Rosen or Greenberger-Horne- 
Zeilinger state comprised of Josephson junctions. It may 
also be feasible to perform a "quantum nondemolition" 
measurement on a one-photon resaaator much as is done 
using a two-level atom in a cavityO We hope to develop 
these ideas in a future publication. 

The remainder of this paper is organized as follows. In 
Section II, we introduce our model Hamiltonian and its 
physical realization; we also outline the method used to 
solve for its eigenstates. Section III contains our princi- 
pal results. We describe the reduced Hamiltonian, and 
we compare the time evolution of the system as predicted 
by the full Hamiltonian and by the two-level approxima- 
tion. In the final section, the main points of the paper 
are summarized, and some of the physical parameters 
entering the calculations are estimated. 

II. MODEL HAMILTONIAN AND ITS 
DIAGONALIZATION 

A. Hamiltonian 

We consider a group of underdamped Josephson junc- 
tions in a lossless electromagnetic cavity which can sup- 
port a single photon mode of angular frequency lu. The 
system is assumed to be described by the following model 
Hamiltonian: 

H = W P hoto n + X) W J (■?')• (!) 

3 

Here 7i p hoton is the Hamiltonian of the cavity mode, 
which we express in the form 7i p hoton = hu) (at a + 1/2) , 
where at and a are the photon creation and annihilation 
operators with commutation relation [a, at] = 1. Ti.j(j) 
is the Hamiltonian of the jth junction, which implicitly 
contains the coupling to the cavity; it can be written as 

Hj(j) = \u {n 3 - n.jf - Jcos 7i , (2) 

where the first and second terms represent the charging 
and Josephson energies of a junction. U = 4e 2 /C is the 
junction charging energy, Ej = hl c / (2e) is the Josephson 
coupling energy, C is the effective junction capacitance, 
and I c the junction critical current, rij is an operator 
which represents the difference between the number of 
Cooper pairs on the two grains comprising the junction, 
and fij is related to the gate voltage applied across the 
junction. Although fij could, in principle, be controlled 
separately for each junction, we assume that all the n^'s 



are equal. The last parameter, jj, is the gauge invariant 
phase of the junction, as defined below. We assume that 
all the junctions have the same J and U, and also that 
J and U are comparable in magnitude. Furthermore, we 
neglect the current through the shunt resistance. Finally, 
we consider only junctions such that both J and U are 
small compared to the superconducting gap A, which is 
thus the largest energy scale in the system. 

The interaction between the junctions and the electro- 
magnetic field is contained in the gauge-invariant phase 
difference 

7i = ^- - (2tt/$o) /A(x,t)-#- (3) 

J 3 

Here <j)j is the phase difference across the junction in a 
particular gauge; it satisfies the commutation relation 
[nk,4>j] = —i5jk- $0 = hc/(2e) is the flux quantum, 
A(x, t) is the vector potential, and the line integral is 
taken across the ]th junction. We consider only a vec- 
tor potential arising from the electromagnetic field of the 
cavity normal mode. We also assume that the junctions 
are all located within a distance small compared to the 
radiation wavelength, so that A(x) is approximately uni- 
form throughout the region of the junctions. If the junc- 
tions all have the same orientation, this condition implies 
that each junction has the same coupling to the cavity 
electromagnetic field. The generalization to nonuniform 
fields is straightforward, however. 

In Gaussian units, and in the Coulomb gauge 
(V • A = 0), the vector potential takes the form 
A = y/hc 2 /u}V( a + aM e, where e is the unit polariza- 
tion vector of the cavity mode, and V is the volume of 
the cavity. Thus we can define the coupling g such that 
(27t/$ ) fA-di= (g/V2)(a + at), or 

a = 4 e A /— ^— e • I (4) 

where £ is the thickness of the insulating layer in 
a junction. With the change of variables p = 
iy/hu/2 (at — a) ; q = ^Jh/{2uj) (at + a) , jj can be 
rewritten 

7j = h ~ 9 VW^q, (5) 
where p and q satisfy the canonical commutation relation 
\p, q] =-ih. 

Finally, we decompose the Hamiltonian ([lj) as follows: 

W = Wphoton + ^ ftjij) + ^int, (6) 
3 

where 

H°j(j) = ±U(n j -fi j ) 2 -J cos fr, (7) 

and 

Hint = -J ^(008 7^ - COS$j). (8) 
3 

In this form, the photon-junction interaction is entirely 
contained in the last term. 
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B. Method of solution 

As in Rcf. |l], we diagonalize H in a complete ba- 
sis formed from the direct product of the eigenfunc- 
tions of Wphoton and TCj(j). The eigenf unctions of 
Hphoton are, of course, harmonic oscillator eigenstates; 
the normalized nth eigenstate has wave function h n (q) = 

, 1 exp(—y 2 /2)H n (y), where H n {y) is a Hermite 

V V^ 2 " n! 

polynomial of order n, and y = (lj/K) 1 ' 2 q. An eigen- 
function ipi&j) °f ^j(j) with eigenvalue E^\fij) can be 
written as ipi^j) = exp(i fij <f)j) rj{(pj), and satisfies the 
Schrodinger equation, H®(j) ip(4>j) = E^' ip((j>j)- Using 
the representation rij — —id/d(j)j, which follows from 
the commutation relation [rij,(j)j\ = —i, we write the 
Schrodinger equation as 

d 2 Y(x) /8E<~ti \ 
^^+(^ r + 2Qcos2x)Y(x) = 0, (9) 

where x = <f>j/2, Y(x) = n(^-/2), and Q = 4J/U. 
This is a Mathieu equation with characteristic value 
a U) = SE^/U, and potential of strength Q. The 
eigenvalues E^ are determined by the requirement that 
^(^i) = i'i'fij + 27r), or equivalently Y(x + n) = 
exp(— 2 i fij 7r) Y{x). Thus we are interested only in the 
Floquet solutions Yk j (x) of the Mathieu equation (||), 
with Floquet exponent Vj = 2 kj — 2 fij , where kj = 
0,±1,±2,.... We denote the corresponding eigenvalue 
of 7Yj(j) by Ef?'(fij). It is sufficient to consider only 

< fij < 0.5, because the eigenvalues of Hj(j) are pe- 
riodic in fij with period unity, and are symmetric about 
fij = 1/2. 

Now, any eigenstate ■■■4>ni Q)i of the Hamilto- 

nian Ti can be expressed as a linear combination of prod- 
uct wave functions consisting of eigenfunctions of 7i p hoton 
and the Hj(j)'s, i. e. , 

N 

$((f>i,...<f>N,q) = ^2 A kl ... kN!n h n (q) Y[ip kj (ct)j), 

. ( 10 ) 

where ipkjifij) = exp(i rij 4>j)Yk 3 \(j)j/2), and N is the 
number of junctions in the resonator. The only term in 
Ti. which is not diagonal in this product basis is the inter- 
action term 7Yi n t - We can obtain the product coefficients 
and eigenvalues by directly diagonalizing the Hamilto- 
nian of Eq. (Q) in this product basis, i. e., by solving 
the Schrodinger equation, H ^((f>, q) — E^((f),q), where 
we use the shorthand notation (f> = (</>i, ...0tv)- Once the 
eigenvalues E and expansion coefficients A kl ,,± Nin are 
known, we can study the time evolution of the system 
given its state at time t = 0, and hence obtain the time- 
dependent expectation value of any operator O at time t 
using 

(0(t)) =Tr [p(t)0]. (II) 

Here p(t) is the density matrix of the system at time 
t, defined by p(t) = U (t) p(0)W (t) , where U(t) = 



exp{—iHt/h) is the evolution operator, and p(0) is de- 
termined by the initial state of the system. 

We have carried out extensive calculations for the 
Hamiltonian (1), using the product basis described 
above. In our calculations we have arbitrarily used 
Q = 0.7 and Ku/U = 0.3, representing a case where the 
charging energy, Josephson coupling, and photon energy 
are all comparable. The attainability of these parameters 
in realistic systems is discussed further below. Most of 
our calculations have also been carried out near the value 
n = 0.258, corresponding to the case where the energy 
difference between the two lowest junction levels equals a 
single quantum of photon energy. In our single-junction 
calculations, El this choice leads to resonant absorption of 
a single photon and to maximum entanglement between 
the junction and the photon states. For larger number 
of junctions ( N > 4) we found that we have to tune n 
slightly away from the resonance condition of the single 
junction to get the system on resonance again. 

In all of our calculations, we use a truncated basis of 
the eigenstates of Hj(j), and of Hphoton- For N = 1, 2, 3 
junctions, we include the first nine photon eigenstates 
and the first five levels of each junction, i. e., a total 
of 9 x 5 W states. For N = 4, we use fewer eigenstates: 
the lowest eight photon states and the first four junction 
levels. In all cases, we confirmed the numerical stabil- 
ity of our results by increasing the dimensionality of the 
basis and finding that there is almost no change in our re- 
sults. As another check, we have calculated Tr(p), which 
should be time-independent and equal to unity for our 
normalization. We find that this trace does remain equal 
to unity to high accuracy for all times considered, thus 
providing evidence that we have correctly diagonalized 
the Hamiltonian. 



III. RESULTS AND DISCUSSION 

A. Reduced Hamiltonian 

The main goal of this paper is to study how well the 
results of this exact diagonalization procedure can be 
matched by a reduced Hamiltonian which includes only 
the lowest two energy levels from each junction. For the 
one-junction case, the two-level approximation was al- 
ready tested in Ref. [l] at h — 0.258. At this value of 
n, the splitting between the two lowest levels equals the 
single-photon energy Tiu). The reduced Hamiltonian then 
takes the Jaynes-Cummings formfl 

H JC = hLj(^a + S z )+t^ S-+£*aS+. (12) 

Here the junction is represented as a quantum spin-1/2, 
with = 2S Z , and the junction energy zero is 

shifted so that the lowest two levels are ±hui/2. The ma- 
trix element coupling the junction and the field is given 
by £ = (/|iWi n t|i}, where Hint is the one-junction ver- 
sion of expression (||); the ket \i) = \n — l;fc = 0) is 
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a wave function containing one photon and the junc- 
tion in the ground state (k — 0), while |/) = \n = 
0; k = 1) is a state with zero photons and the junction 
in its first excited state (k = 1). To first order in g, 
Hint ~ —Jg-\/u}/hqsm(j), and hence, to the same order 
in g, 

-(Jg/V8) (k 



(a) 



Olsin^l/c = 1) 



(13) 



This completes the definition of the reduced Hamiltonian 
for one junction, and establishes its relation to the full 
Hamiltonian of Eq. (|l]) . 

A naive extension of the reduced Hamiltonian ( |l2| ) to 
N junctions would be 



tl D = tin + i^T 



where 



A' 



(14) 



o/O 

tl n 



N 



Tin = 



Here S_ , SV - ', and arc the lowering, raising, and 
inversion operators of the jth junction, which satisfy the 
commutation relations 



rc(j) Q W 



± 



[5^,5i fe) ] 



The Hamiltonian (Q) resembles the one-mode Dicke 
model,B£3 which describes N identical two-level atoms 
within a radiation wavelength of each other, all inter- 
acting with a single mode resonator. The initial state 
problem considered later in this section has been solved 
exactly for this Hamiltonian in a cavity conta ining one 
or more accessible electromagnetic modes.E3E3 

For the present system, when the resonator con- 
tains more than one photon, the Hamiltonian (jl4|) docs 
not reproduce the full time-dependent solution to the 
Schrodinger equation for N > 2, at least for our nu- 
merical parameters. Nevertheless, as shown below, the 
full numerical solution is still well approximated by an- 
other reduced Hamiltonian in which each junction is rep- 
resented by a basis of two levels only, namely 



tlD' — tln> + tin', 



(15) 



where the zeroth order Hamiltonian H^ D , is the same as 
the Dicke TiP D appearing in Eq. (|l4|), but the interaction 
part of the Hamiltonian Tip, is now 



N 



0) 



3 = 1 



N 



nJ2 (sV>s!?> + s®s w 



(16) 



7,fc = l 
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FIG. 1: (a) Time-dependence of (H°) (full curve) and of 
(^photon) (dashed curve), plotted for one junction with n = 
0.258 (on-resonance condition). The system is prepared at 
time t = in the state |n = 1 ; k = 0). The other parameters 
are Q = 0.7, cu/U — 0.3, and g = 0.15. (b) Same as (a) except 
that we show the time-dependent probability for finding the 
Josephson junction (full curve) and the photonic resonator 
(dashed curve) in their first excited states. 



The physics behind this reduced Hamiltonian0 is dis- 
cussed further below, in Sec. III. D. 2. 



B. Review of single junction 

We first review the time-dependent solution to the. 
Schrodinger equation for a single junction in a resonator .El 
If the initial junction-cavity state is \n = 1 ; k = 0), then 
the system evolves in time according to the evolution op- 
erator U(t) = exp(— iTit/fi). Fig. 1 shows the resulting 
time-dependent expectation values (Tij) and (7i p hoton), 
as calculated using Eq. (|ll|) and the full Hamiltonian 
H. Fig. 1 also shows the time-dependent probabilities 
for finding the junction and the resonator in their first 
excited states, given this initial state. The system wave 
function is evidently an entangled state oscillating be- 
tween \n = 1; k = 0) and |n = 0; k = 1), with a pe- 
riod T w 675H/U for our particular choice of g, Q/U 
and huj/U. When the Rabi period for oscillations be- 
tween the Josephson-photon entangled states is calcu- 
lated using the two-level approximation, it is found to be 
7r?i/|£| = 667 fb/U, in almost perfect agreement with the 
result obtained from the full Hamiltonian. 



C. Several junctions, asymmetric initial state 

Next, we consider N junctions such that initially (at 
time t = 0) one junction is in the first excited state, while 
the remaining (N — 1) junctions are in the ground state, 
and the cavity initially contains zero photons. Thus all 
the junctions interact with the vacuum fluctuations of 
the cavity mode via the Hamiltonian Ti. m t . 

Our primary interest is to compare the time evolu- 
tion calculated from the exact Hamiltonian with that 
obtained from the two reduced Hamiltonians discussed 
above. We will focus our attention on the "inversion" 



5 



1 

0.9 
0.8 
0.7 
0.6 
0.5 



(a) N=2 




0.5 



1.5 



1 

0.95 
0.9 

0.85 
0.8 

0.75 



(b) N=4 




0.5 



1.5 



Time (10 fi/U) 



FIG. 2: The time-dependent inversion S(t) [Eq. (17)] for 
(a) N — 2 and (b) N — 4 junctions in a resonator, calculated 
using the full Hamiltonian. Initially, one junction is in its first 
excited state, while all the other junctions and the resonator 
are in their ground states. 



S(t) of the system, defined as 
S(t) = (S z (t) 



(17) 



where the average is taken with respect to the quantum 
state of the system. S(t) is the sum of the probabilities 
that the junctions are in their first excited states, and sat- 
isfies < S(t) < N. This quantity is of interest because 
1iujj§ — I, is the rate at which energy is transferred from 
the junctions to the cavity or vice versa. 

For the initial state considered in this section, it has 
been proventa for the Dicke model [Eq. (14)] that 



S(t)=l-- a m\\Z\VNt/h) 



(18) 



Thus, the frequency of energy transfer between the pho- 
ton field and the junctions is proportional to \/N. The 
group of N junctions behaves somewhat like a single junc- 
tion with a coupling to the cavity mode of y/N £ instead 
of £ . It also follows from Eq. ([l8]) that the photon field 
never gains more than 1/N of the energy initially stored 
in the excited junction. Thus, as TV increases, it be- 
comes progressively less likely that a photon will be emit- 
ted into the cavity, the energy remaining largely trapped 
within the junctions ( actually, within the initially ex- 
cited junction). The existence of N — 1 unexcited junc- 
tions (or "atoms" ) has the effect of preventing the emis- 
sion of a photon from the initially excited juncjtipn,E3 a 
phenomenon known as "radiation suppression. "123 

When we calculate S(t) using our direct diagonaliza- 
tion scheme, we find that the numerical results are in 
excellent agreement with the predictions of the two-level 
approximation ( |l8| ) for N = 1,2,3 and 4 junctions in the 
resonator, at least for several Rabi periods of oscillation. 
In Fig. 2, we show the inversion as calculated using the 
full basis for N = 2 and N = 4. 



D. Dicke symmetric states 



which the initial state of the system is symmetric under 
the interchange of any two junctions. We have solved 
for the time-dependent state of the system using the full 
Hamiltonian, and compared the results to the approxi- 
mate Hamiltonians ( |l4| ) and (|l5|). Results of our calcu- 
lations are shown in Figs. 3-5. 

As can be seen from these figures, the one-mode Dicke 
model defined in Eq. ( |l4| ) does not reproduce the full 
numerical solution for n > 2. But the numerical solu- 
tion is well approximated by the reduced model Hamil- 
tonian ([To]), in which each junction is still represented 
by only two levels. This Hamiltonian differs from the 
Dicke model only in the last term, which is proportional 
to f2 and which represents an effective direct interaction 
between the junctions. This term is needed in order to 
obtain a good fit between the predictions of the reduced 
model and the exact Hamiltonian, even though the junc- 
tions in the original model Hamiltonian [Eq. ([!])] have 
only an indirect coupling with each other through the 
electromagnetic field. 

A similar model Hamiltonian to Eq. ( |l5| ) has been 
studied in Ref. [2^ for two atoms in a cavity. This model 
differs from the present one, however, in that the di- 
rect coupling in Ref. 22 represents a specific term in the 
Hamiltonian, namely, an effective dipole-dipole interac- 
tion between the atoms. In the calculations to be pre- 
sented below, the coupling strength Q in the Hamiltonian 
of Eq. ([l5]) is treated simply as a parameter determined 
by best fitting to the numerical data. However, it can 
also be derived explicitly from the full Hamiltonian, as 
we show below. 

There is one further point to be made about these 
numerical results. Namely, so long as n < N, we find 
that the interaction term £ which best fits the computed 
time evolution is numerically equal to that which enters 
the single-junction reduced Hamiltonian [Eq. (|l3|)] when 
there is no more than one photon present. This result is 
intuitively reasonable, because in this symmetric case all 
the junctions in the resonator share equally in the pho- 
tons, and thus, for n < N, none will absorb more than 
one photon. 

Now we will discuss the basis used to study the Hamil- 
tonian (|l5|). The Hilbert space of N two- level junctions 
is 2 N dimensional. Thus the system is analytically in- 
tractable even for a small value of N unless there is some 
symmetry of the Hamiltonian which confines the evolu- 
tion of the system to some lower invariant subspace. The 
Hamiltonian ( |l5"| ) has two useful symmetries of this kind. 
First, a symmetric initial state is an eigenstate of the per- 
mutation operator Vjk, which exchanges the jth and fcth 
junctions. Since Vjk commutes with the Hamiltonian of 
Eq. (|l5|), a symmetric initial state evolving under H.D' 
will remain symmetric at later times. The second sym- 
metry involves the excitation number operator 



Next, we consider an initial state in which all the junc- 
tions are in the ground state and there are n photons 
in the system (n — 1, 2, 3). This is a simple case in 



N 

M=a^a + Y / i 



(19) 
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FIG. 3: (a) Time-dependent probability for finding any one of the junctions (full curve) or the photonic resonator (dashed 
curve) in its first excited state for N — 2 (top) and N — 4 (bottom). The system is prepared at time t = such that there is 
one photon in the system, and all the junctions are in their ground states, (b) Same as (a) except we show the inversion S(t) 
of the system as defined in Eq. (fLTl) . 



which also commutes with 7irj', and thus is also a con- 
stant of the motion. In this case it is convenient to study 
the Hamiltonian TLd' using the basis vectors formed from 
the direct product of the eigenstates of M. and a^a. We 
denote these states as \m:n), where 



A4\m;n) — m\m;n), 
a*a\m\n) = n\nv,ri), 



(20) 



and < n < m. Note that the time evolution of the 
system is restricted to a (m + l)-dimensional .subspace 
spanned by the basis vectors |m; 0), . . . , |to; to) £3 

It is convenient to write out the state \m;n) explic- 
itly. There are n photon excitations in the cavity, and 
to — n excitations distributed symmetrically among the 
N junctions. Thus we may write \m; n) = \n) ®\m — n) j, 
where \n) is a Fock state and |m — n)j is a normalized 
symmetric Dicke state of the junctions, described by 



I*). 



k\(N-k)\ 



-1/2 



eih+>* n i->*». 

V 1=1 m ^l...fe 

( 21 ) 

where = \h = 1) (|— )j, = \h = 0) ) is the spin up 
(down) state of the /th junction, and the summation is 
over permutations V among the N junctions. 

We now discuss the cases of m = 1,2,3 excitations, 
and compare the time-dependent states obtained from 
the full numerical solution to the predictions of the two 
level approximation, using the basis \m;n). 



1. One excitation (m — 1) 

In this case, our numerical results are well described 
by the Dicke Hamiltonian [Eq. (jlj)], at least for N < 
4. This behavior is expected, since the junctions cannot 
interact directly with each other in the presence of only 
one photon. For one excitation in the system, an exact 
solution for a twp-level system can be found even in a 
damped cavity.Cj In the absence of damping, the wave 
function of the system in the interaction picture can be 
written in the two-level approximation as: 

\V\t)) = C 1 (t)\l;l) + C (t)\l;0). (22) 

From the time-dependent Schrodinger equation (taking 
fi = 1 henceforth), 



we obtain 



K\^ I (t)) = i -\^(t)), 



(23) 



dt 
.dC (t) 

l : 

dt 



C VNd(t). 



Taking the junctions initially to be in their ground states 
and assuming n = 1 photons in the system (i. e. 
Ci(t)\t=o = 1), we find after some algebra that 

Ci(t) 
C (t) 



cos(|£| VNt), 
sin(|£| VNt). 
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FIG. 4: Time-dependent inversion S(t) for m — 2 excitations and TV = 2, 3 and 4 junctions in the resonator [(a), (b) and (c) 
respectively] . The initial state has two photons in the resonator, and all the junctions in their ground states. The solid line is 
the full numerical solution, and the dotted line follows from Eqs. (j^j. Note in particular the decrease in the Rabi period of 
oscillation as TV increases. Other parameters are the same as in Fig. 1 except that fij — 0.2598 for TV = 4. The Dicke model 
( not shown) would give a sinusoidal oscillation with no beats, both here and in Fig. 5. 



Thus the system oscillates between the states \m = 1; n = 
1) and \m = l;n = 0)- that is the junctions alternately 
absorb the photon of the field, then re-emit it, in a re- 
versible evolution, as expected in a lossless cavity. The 
oscillation period T oc TV -1 / 2 . From this solution, the 
inversion of the system [Eq. (|l^)] is 



(24) 



S(t)=faa 4 (\Z\VNt), 
and hence the rate of radiation into the cavity is 
dS(t) 



lit) 



dt 



Lo\£\VNsm(2\Z\VNt). (25) 



Both the oscillation frequency and the amplitude of the 
rate of radiation are proportional to \/~N, a signature of a 
cooperative phenomenon. For very short times, the rate 
of radiation into the cavity is 



I(t) sa -2uj\£\ 2 Nt 



(26) 



- that is, TV times the rate of a single junction. Thus 
the presence of TV-j unctions in the cavity gives rise to a 
collective behavior. 

Fig. 3(a) shows the time-dependent probability that 
the Josephson junctions (full curve) and the photonic res- 
onator (dashed curve) are in their first excited states for 
TV = 2 and TV = 4. The system is prepared at time t = 
with all the junctions in the ground state, and the res- 
onator in the number state \n = 1). The system develops 



into an entangled state such that the photon is shared 
equally among all the junctions. This behavior can be 
seen both in the two-level approximation and in the re- 
sults of an exact calculation. In Fig. 3(b) we show the 
inversion, S(t), calculated using the full Hamiltonian for 
N = 2 and TV = 4 junctions and one excitation (m = 1). 
The numerical results are clearly in excellent agreement 
with the predictions of the two-level approximation [Eq. 

2. Two excitations (m — 2) 

For m = 2 we need consider only three eigenstates of 
the Hamiltonian. In the interaction picture the system 
wave function is 

\*\t)) =C 2 (t)\2; 2) + d(t)|2;l) + C (i)|2;0), (27) 

where the coefficients satisfy 

. dC 2 (t) 



dt 
dCx{t) 

dt 
dC (t) 

dt 



= C V2NC 2 (t) + n d(t) +£y/2(N-l) C (t), 



r y/2(N-l)d(t). 



Assuming an initial state |2;2), i. e. a state containing 
two photons, with all the junctions in their ground states, 
we find that the solution of these coupled equations is 



C 2 (t) = c{2N _ i) {aN-l) + Ne-^[C cos(C t) + i n sin(C t)} } , 

Cx{t) = -iy/2N(Z/Q e-? int sm((t), (28) 
Co{t) = 2(2N-l)C I 6 '" ^ [2C C ° S(a) + ^ Sln(a)] ~ 2C I ' 
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where C — 5 y/ ^ 2 + 8 \^\ 2 (2 N — 1). The probability of a two-junction excitation of the system is P2(t) 
which can be reduced to 

N(N - 1) 



(2N- 1) 2 C 2 



2C 2 - 2 (2 JV — 1) |£| 2 sin 2 (Ct) 

- 2 C 2 cos(C t) cos(n t/2)-Q( sin(C i) sin(fi t/2) 
I 



|Co(t)| 2 , 



(29) 



For fl = 0, this becomes 

P2{t) = {2N-\f t 008 ^)- 1 ] 2 ' ( 3 °) 

so that the Rabi period of oscillation is now cu 2 — [2(2 N— 
I)] 1 ' 2 Since \P 2 (t)\ < N(N- 1)/(2N- l) 2 < 1/4, the 
system is unlikely to contain two excitations simultane- 
ously. Even if fl ^ 0, the same statement is still valid, but 
the time evolution is no longer perfectly sinusoidal, but 
instead is characterized by beats with a period inversely 
proportional to tt. 

The inversion of the system is given by 

S(t) = \C 1 (t)\ 2 +2\Co(t)\ 2 . (31) 

A full expression for the inversion readily calculated from 
Eqs. (H|) is quite complicated, but is not given here. Fig. 
4 shows the results obtained numerically for N — 2,3, 
and 4, using the full basis (solid line), as well as the 
inversion calculated from Eq. (|3l]) using the modified 
Dicke model with £1 ^ (dashed line). 

We now discuss the origin of the extra direct interac- 
tion in the effective Hamiltonian (|ll|), in light of these 
numerical results. Eq. (^) expresses the wave function 
for two excitations as a linear combination of three states, 
which we will denote compactly as |a), \b), and \c) respec- 
tively. In the absence of Hint, these states are degenerate. 
According to the Dicke model, Hint breaks this degener- 
acy because some of the matrix elements of Hint between 
these states are nonzero. This is the first- order effect of 
Hmt (or, equivalently, in g). 

However, there is also a second-order effect. Namely, 
as is well known frcwa degenerate second-order Lowdin 
perturbation theorytJ the effective Hamiltonian matrix 
elements between, for example, states (a\ and \b) should 



be written 

(a\H eS \b) = (a\Hi nt \b) + ± (^f)(f^) _ (32) 

r tja ~ hr 

Here the sum runs over all states \r) other than \a), \b), 
and I c) . These are the terms which give rise to the extra 
interaction proportional to in Eq. ( |l5| ) . 

To confirm this hypothesis, we have calculated these 
matrix elements, including the terms second order in 
Hint, through order g 2 . We do find that these terms 
produce a direct interaction between the junctions of the 
form (|l5|), as found numerically. Specifically, the matrix 
element (b\H e s\b) becomes nonzero when these second- 
order terms are included. Furthermore, we find that the 
value of f2 found from this direct calculation is com- 
parable to that found numerically. Specifically, for the 
the cases mentioned above, we find VL/{\£\ 2 /lu) — 4.89, 
9.65, 14.42 for N — 2, 3, and 4. For comparison, the 
best fit values of f2 to the numerical results of Fig. 4 are 
fl/(\^\ 2 /Lj) w 5.74,4.46, 13.40 for N = 2, 3 and 4 respec- 
tively. 



3. Three excitations (m = 3) 

In this case the wave function at any time t in the 
interaction picture is 

|* I (t))=^C i (i)|3;i), (33) 

where the expansion coefficients satisfy 



dC 3 (t) 

dt 
dC 2 (t) 

dt 
dd(t) 

dt 
dC (t) 
dt 



ZV3NC 2 (t), 



C V3N C 3 (t) + nC 2 {t) + 2ZVN-lC 1 (t) , 



2? VN~lC 2 (t) + nC\{t)+£ ^3 (JV-2)Co(t), 
r V3(JV-2)Ci(i). 



The inversion of the system in this case is 

S(t) = \C 2 (t)\ 2 + 2\C 1 (t)\ 2 + 3\C (t)\ 2 . 



(34) 



r 



In Fig. 5 we show the results obtained from the full nu- 
merical solution (full line) and from Eq. (p4|) (dashed 
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(a) N=3 



(b) N=4 



(a) N=1 



(b) N=2 
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FIG. 5: Same as Fig. 3, except that the initial state has three 
photons in the resonator, and calculations are carried out for 
N = 3 and N = 4 junctions, as indicated in the figure. 

line) for N = 3, and N — 4 junctions. Here the best 
fit values of are given such that Sl/(|^| 2 /a;) « 7 and 
13.4 respectively. For comparison, the second-order de- 
generate perturbation calculation mentioned above gives 
fl/(\£\ 2 /uj) = 10 and 15.5. 

E. Coherent state 

As a last non-trivial case we consider the system such 
the resonator is initially in a coherent state. This case is 
of special interest because the number eigenstates of the 
electromagnetic field are generally difficult to generate. 
Nevertheless, as we will see below, our main conclusion, 
namely the enhancement of the junction-resonator inter- 
action with increasing N still holds in this case. 

A coherent state \f3) with an average number of pho- 
tons (n) = \f3\ 2 can be written as a linear superposition 
of number states \s) with a Poisson distribution: 

|^) = ^ e -|^/ 2 ^l| s) . (35) 

For N = 1, it is easily proved that the inversion of the 
system is given by: 

5(t) = i-ie-l^ 2 f;^co S (2e^) ! (36) 

which shows the phenomenon of " collapse and revival" 
i. e. the envelope of the Rabi oscillations periodically 
"collapses" to a half with a subsequent "revival" at later 
times£n23 For N > 1, the analysis of the previous subsec- 
tion can be readily generalized to treat a coherent initial 
state. 

As an example, we show in Fig. 6 the inversion S(t) 
for N — 1 and N — 2 assuming (3 — \/2. The solid 
lines follow from the full Hamiltonian while the dashed 
lines are from the reduced Hamiltonians of Eqs. ( p^ ) 
and (|la). The best fit value of f2 in this case is ft « 
5.74 This value of f2 is the same as for the number- 

state considered before with m — 2 excitations. The 
enhancement of the junction-cavity interaction can be 
seen from the decrease of the period of oscillation in Fig. 
6(b) compared to Fig. 6(a). 




012345 012345 



Time ( 10 3 hi U) 

FIG. 6: The time-dependent inversion S(t) for (a) N = 1 
and (b) N = 2 junctions given that the initial state of the 
resonator is a coherent state with average photon number 
(n) = 2, and the junctions are initially in their ground states. 
The solid line denotes the full numerical solution, and the 
dotted line follows from the two-level approximation. 



IV. SUMMARY AND CONCLUSIONS 

The present paper extends our previous publication!]] to 
the case of N > 1 identical junctions interacting with a 
single lossless cavity mode. Although oversimplified, this 
model has the advantage that we can obtain detailed, 
numerically exact dynamical solutions with a range of 
initial conditions. Our treatment of the Hamiltonian is 
entirely quantum-mechanical. 

We have diagonalized the Hamiltonian exactly, using 
a basis large enough to insure accuracy, and compared 
our results to those of a reduced Hamiltonian [Eqs. ( |l4| ) 
or ( [llf ) ], in which each junction is represented by two 
levels only. The reduced Hamiltonian resembles the one- 
mode Dicke Hamiltonian but contains a direct interac- 
tion term between the junctions. For short times (of the 
order of several Rabi periods) the agreement between re- 
sults obtained from the two-level Hamiltonian of Eqs. 
dl4| ) or (JlJ) and those found from the full Hamiltonian 
[Eq. ([!])] is excellent, but decreases somewhat for longer 
times. Also, as expected for a Dicke-like model, the sys- 
tem of several identical junctions behaves just like a single 
junction for short times, but with an enhanced junction- 
cavity coupling strength £ y/N, and a correspondingly 
decreased Rabi period. This reduced Rabi period might 
be smaller than the coherence time, making these oscil- 
lations experimentally observable. 

The most striking feature of our numerical results is 
the presence of an effective direct interaction between 
the junctions resulting from the cavity. This interaction 
has the form of a dipole-dipole interaction between the 
"spins" representing the two lowest levels of each junc- 
tion. This extra interaction must be included in the effec- 
tive two-level Hamiltonian in order for it to reproduce the 
numerical results obtained from the full Hamiltonian at 
long times. While the analogous term in atomic physics 
is a true dipole-dipole interaction, in our model it arises 
through the effects of the higher energy levels on the 
degenerates states, as calculated through second order 
Lowdin degenerate perturbation theory. 

Before concluding, we briefly discuss the case of non- 
identical junctions. Obviously, this is an important con- 
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sideration for Josephson junctions since, in contrast to 
atoms, it is impossible to fabricate perfectly identical 
Josephson junctions. To test whether the coherence 
survives when the junctions are non-identical, we have 
repeated the two-junction calculations shown in Figs. 
2-4 using junctions which differ slightly from one an- 
other. For convenience, we considered two junctions 
with slightly different values of J, but the same values of 
U. With these parameters, the resonant values of n are 
slightly different for the two junctions in the same cav- 
ity. We find that, if the two J's differ only slightly (by 
only 1-3% for our choice of J/U), the coherent behavior, 
as signaled by the y/2 enhancement of the Rabi oscilla- 
tion frequency, survives. For greater differences between 
the J's, the coherence disappears and one finds instead 
two distinct resonances, each with an unenhanced Rabi 
oscillation frequency. 

We have also carried out similar calculations in other 
parameter regimes - specifically, for Q varying between 
0.01 and about 3 (with g remaining unchanged). The 
qualitative results - that is, the presence of Rabi oscilla- 
tions and the characteristic frequency changes associated 
with coherence - are unchanged. 

Finally, we briefly discuss whether the present system 
could be constructed in the laboratory. The values of 
the parameters entering the model have been estimated 
previously!!! For a cylindrical cavity of length d and radius 
r, a typical g [ Eq. (||)] was found to be g ~ 3£ y/a/(rd), 
where a = e 2 /(hc) ~ 1/137 is the fine structure constant, 
and £\\ is the junction thickness. Achievable values of the 
other parameters were estimated as U ~ 10 -15 erg ~ 8K, 
and u) ~ 2c/r. In order to attain Q ~ 0.7, as used in our 
calculations, we need J/Ub ~ 3 K, or a critical current 



I c ~ 0.1/iA. To obtain huo/U ~ 0.3, these estimates 
give x ~ 10 -3 cm and r ~ 0.1—1 cm, where a; is a 
typical junction dimension (e. g., island size, junction 
radius). These values, which would give uj ~ 100 GHz, 
are achievable in present technology. The cor resp onding 
maximum power radiated into the cavity [Eq. |25[ will be 
of the order of lO" 11 ^ W. 

However, the value g — 0.15, used in our calculations, 
is much larger than the value g ~ 10~ 5 which was previ- 
ously estimated 1 ! to result from the parameters just de- 
scribed. Such a small g would certainly cause the two- 
level approximation to be even better than found in our 
calculations here. But it would have the disadvantage 
of producing longer Rabi oscillation time, and hence, of 
making coherence over a Rabi period more difficult to 
achieve; it would also lead to a much lower power radi- 
ated. Thus, the largest possible g is clearly desirable in 
order to detect this effect. The Rabi period could also 
be decreased by junctions with larger critical currents 
than those considered here. However, the critical current 
is limited by the requirement huj <C A, where A is the 
superconducting gap. Thus, the best hope for observ- 
ing these oscillations may be to reduce the decoherence 
effects. 
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